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I. INTRODUCTION 



Grand unification [1] of tlie standard model gauge group into a simple group with the 
matter and Higgs particles transforming under irreducible representations of this simple 
group is an attractive framework for pursuing physics beyond standard model (SM). On the 
other hand supersymmetry [2] is at present the only framework in which the Higgs sector 
of the Standard Model (SM), so crucial for its internal consistency, is natural. A much 
favored implementation of the idea of low energy supersymmetry (SUSY) is the Minimal 
Supersymmetric Standard Model (MSSM) [3] , which is obtained by doubhng the number of 
states of the SM, and introducing a second Higgs doublet (with the opposite hypercharge 
of the SM Higgs doublet) to generate masses for all the fermions and to cancel triangle 
gauge anomalies. The MSSM leads to a successful prediction [4] for the ratios of the three 
gauge couplings, with a unification scale of Mq — 10^^ GeV and supersymmetric thresh- 
olds near the electroweak scale (Mg ~ 1 TeV). The possibihty of such a unification scale 
in the supersymmetric context would naturally suggest the framework of supersymmetric 
grand unification at the scale Mq- Furthermore, the presence of supersymmetry would also 
make the hierarchy between the disparate scales associated with the grand unification (Mg), 
and the weak scale Mw ~ 100 GeV technically natural. However, despite many years of 
dedicated experimental work, no direct evidence has been found to vindicate the idea of 
unification, especially its most heralded prediction of proton decay. All the present experi- 
mental data is, however, consistent with the idea of low energy supersymmetry with a grand 
unified scale being around 10^^ GeV. It is therefore important to keep an open mind toward 
all possibilities, conventional and unconventional, in our search for unification of strong and 
electroweak interactions, since there is going to be a concerted effort to search for collider 
signatures of such possibilities at the next generation of experiments [5] . 

At the theoretical level, the gauge group SU{3) x SU{2) x U{1) remains a completely 
unexplained feature of the standard model. However, it can be elegantly unified into simple 
groups like S'f/(5) [6] or 5*0(10) [7]. The 5*0(10) group is particularly appealing, since it has 
a 16 dimensional representation which is large enough to accomodate an entire generation 
of standard model fermions, and can further accomodate a right-handed neutrino as well. 
The 5*0(10) framework can readily be extended to include supersymmetry. Furthermore, a 
complex 10 dimensional representation of 50(10) can be employed to accomodate the two 
Higgs doublets of the low energy minimal supersymmetric standard model. 

Unlike the case of 5*^7(5), the rank of the group 5*0(10) exceeds that of the standard 
model by 1. The breaking of the additional U{1) factor in 5*0(10) group leads to char- 
acteristic non-universal contributions to the otherwise universal soft SUSY breaking scalar 
masses of the minimal supersymmetric models. These additional corrections to the soft 
scalar masses are usually referred to as D— term contributions, and lead to non- universality 
for the soft scalar masses . In general D— term contributions to the SUSY breaking soft 
scalar masses arise whenever a gauge symmetry is spontanesously broken with a reduction 
of rank [8]. These D— term contributions have important phenomenological consequences 
at low energies as they allow one to reach certain regions of parameter space which are 
not otherwise accessible with universal boundary conditions [9-11]. Non-universality, and 
in particular the £>— term contributions, may have a dramatic impact on the sum rules [12] 
satisfied by the squark and slepton masses. Such effects are likely to help distinguish between 
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different scenarios for breaking of grand unified symmetry at high energies [13-15]. 

In this work we systematically consider the D— term non-universality that is generated 
in 5*0(10) unification. This is the simplest of grand unification models where a D— term 
non-universality is inevitable. Indeed, since 5'O(10) contains SU{5) x C/(l) as one its "nat- 
ural" subgroups, one common route for breaking SO{10) is via SO{10) SU{5) x U{l)z, 
with SU{5) subsequently breaking via SU{5) SU{3)c x SU{2)l x U{1)x- How- 
ever, in the breaking of 5*0(10) via 5f/(5), there exist two different ways of embedding 
SU{3)c X SU{2)l X U{1)y into SU{5) x U{1). In the "conventional" embedding [7] via 
SU{5), the hypercharge generator Y is identified with the generator X of U{l)x- On the 
other hand, in the "fiipped" [16-18] embedding the hypercharge generator is identified with 
a linear combination of the generators X and Z. We shall study in detail the consequences of 
these two different embeddings of the SM gauge group in 50(10). We should mention here 
that apart from the "natural" subgroup SU{5) x ^7(1), the group 50(10) also has "natural" 
subgroup 50(6) x 50(4). Since 50(6) is isomorphic to 5[/(4) and 50(4) is isomorphic to 
SU{2) X SU{2), 50(10) contains the group 5C/(4) x SU{2) x SU{2). We shall not consider 
this pattern of breaking of 50(10) in this paper. 

The plan of this paper is as follows. In Sec. II, we consider in detail the the term 
non-universality in supersymmetric 50(10) unification. We shall review the embedding of 
the standard model gauge group in the 50(10) group. We shall identify the hypercharge 
generator Y for conventional and fiipped embeddings, and the orthogonal generator Y^ 
whose eigenvalues determines the D— term contributions at the unification scale. In Sec. 
Ill, we write down the renormalization group equations for the sfermion masses and their 
solutions, and derive sum rules for these masses in generality for the embeddings of the SM 
gauge group in the 50(10) group. For the conventional embedding of 50(10), it turns out 
that our sum rules are same as those that were discovered for SU{5) unification with non- 
universal masses [15]. We trace the origin of this to the fact that the D-term non- universality 
in the conventional embedding of 50(10) arises in such a manner that it is degenerate 
within the 10 and 5* representations of SU{5). On the other hand we find the result that 
the sparticle sum rules for flipped embedding, where the D— term non-universality is not 
expected to have a simple pattern, are identical to the case of conventional embedding. By 
considering the nature of the embedding of the conventional and fiipped scenarios is some 
detail we establish why this is the case. A further surprising aspect we find is that the 
individual masses entering the sum rules also do not depend on the manner in which the 
SM is embedded into 50(10). 

In Sec. IV, we present results for the sparticle spectrum based on numerical studies 
and our conclusions. We note that since D term contributions to the scalar masses at 
the unification scale [10] has received considerable attention, our results show that it is 
immaterial whether the embedding is the conventional or the fiipped one. One input we 
have employed is the requirement that the fine tuning problem inherent in the embedding of 
MSSM in 50(10) [19] be alleviated, in both the conventional and fiipped scenarios, by the D 
term contribution, thereby enabling us to fix the sign of this contribution. To our knowledge, 
this is the first time that this has been pointed out in the context of flipped uniflcation. 
Noting that the flne tuning problem is further alleviated by the effects of the right-handed 
neutrino Yukawa coupling to the renormalization group flow from the uniflcation scale Mq 
to a scale Mji ~ lO-*^^-^ GeV, a scale inspired by solar neutrino physics [20], it is possible to 
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render supersymmetric 5*0 (10) unification into a technically natural and viable candidate 
for a complete theory, with unique signatures. 



II. CONVENTIONAL AND FLIPPED EMBEDDING OF THE SM IN SO{W) 

We begin by noting a crucial difference between the SO{10) GUT model and the SU{5) 
GUT model. Whereas the embedding of the SM gauge group into a unified gauge group 
is unique for SU{5), there are two different ways in which the hypercharge group can be 
embedded in 5*0(10). These different embeddings may lead to different predictions for the 
low energy physics. 

In 5*0(10) grand unification, all the matter particles of one family of the Standard 

Model (SM) together with a right handed neutrino belong to the spinor representation 
16. Each such spinor representation 16 can be decomposed under the maximal subgroup 
50(10) D SU{5) X U{l)z as 

16 = 5*3 + 10_i + l_5, (1) 
10 = 52 + 5*_2. (2) 

Furthermore under SU{5) D SU{3)c x SU{2)l x f/(l)x, we have the decomposition 

5 = (3,1)_2 + (1,2)3, 
5* = (3*,1)2 + (1,2)_3, 
10=(3,2)i + (3*,l)_4 + (l,l)6, 

l = (l,l)o. (3) 

We note that C/(l)x, which is the subgroup of 5'C/(5), is not identical with the U{1)y of the 
SM at this stage. Note also that each 16 includes two pairs of (3*, 1) and (1, 1). 

In order to identify the hypercharge group, we consider the decomposition 50(10) D 
SU{h) X U{l)z D 5f/(3)c X SU{2)l x U{1)x x U{1)z. Therefore, the hypercharge U{1)y 
must be a linear combination of U{l)x and U{l)z, i.e. U{1)y C U{1)x x U{1)z- Thus, 
there are two ways to define the hypercharge generator of the SM: 

1. Y = X, 

2. Y^-l(X + 6Z), 

upto an overall normahzation factor. The first case corresponds to the the Georgi-Glashow 
model [6], whereas the second identification of hypercharge corresponds to the fiipped 

case [16,17]. 

In the first case the f/(l) generator of 50(10) that is orthogonal to Y and the diagonal 
generators of SU{3)c and SU{2)l is 

Y^ = -Z, (4) 
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whereas in the flipped case we have for the orthogonal generator^ 

= (5) 

We wiU comment on the normahzation in the following. After a suitable identiflcation of 
the flelds lying in the relevant representations of 5'O(10), the effect of 5'O(10) breaking at 
the uniflcation scale leads to D-term non-universality, which is computed in terms of the 
eigenvalues of the operator Y-^ on the fields which will discuss in the next section. 



III. SOLUTIONS OF RENORMALIZATION GROUP EQUATIONS AND SUM 

RULES 

For the squarks and sleptons of the first and second family (light generations) , the renor- 
malization group (RG) equations for the soft scalar masses are given by 



dt 






(6) 


167r' 

dt 




- i*?^' - P''- 


(7) 


dm\ 
dt 


= 




(8) 


dm\ 
dt 


= -6glM^ - 


■ -^QiMi - -QiS, 


(9) 


167r^ 

dt 


24 

= -f sjAff 




(10) 



where t = Itl^Q/Qq), with Qo being some initial large scale; Ms^2,i are the running gaugino 
masses, ^'3,2,1 are the usual gauge couplings associated with the SM gauge group, and 

S ^ Tr{Ym') = - <^ + E - 2m?^ + - + m?J . (11) 

families 

The U{1)y gauge coupling gi (and cci) is taken to be in a GUT normalization throughout 
this paper. The quantity S evolves according to 

dS _66ai 

~dt - Y4^^ ^^^^ 



which has the solution 



S{t)^S{tG)^y (13) 



In order to determine this, one must work with a set of normalized generators. The expUcit 
normalization factors are adapted from the expressions given in ref. [17]. 
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We note that if yS" = at the initial scale, which would be the case if all the soft sfermion 
and Higgs masses are same, then the RG evolution will maintain it to be zero at all scales. 
The solution for the renormalization group equations (6)-(10) can then be written as 



'"'UL 


= (tc) + C's + C2 + ^Ci + - ^ sin^ ew)Ml cos(2/3) + ^ 




(14) 




1 1 1 

= m^^[tG) +C3 + C2 + —Ci + [-- + -sm UwjMz 


cos(^zpj + 


5 ' 


fl5) 




= ml^itc) + + ^Ci + ^ sin^ 9wMl cos(2/?) - ^K, 






[i-b) 




= {tG) + Cs + ^Ci-^ sin^ OwM^ cos(2/3) + '^K, 






(17) 




- mlfia) + C2 + |Ci + + sin^ ^t^)M| cos(2/?) - 






(18) 




= ml^ito) + C2 + \c, + 008(2/3) - ^X, 






(19) 




= ^L(^g) + Ci - sin^ ^^^Mf cos(2/3) + ^K, 








(20) 



where Ci, C2 and C3 are given by 

Ci(t) = ^ 5.(^)^M,(^)^ i = 1, 2, 3 (21) 



3 3 4 
01 = ^,02 = 1,03=2, (22) 



and 



^ _ q;i(^ 



Q;i(i) 



(23) 



is the contribution of the non-universality parameter 5" to the sfermion masses, and bi — 
-33/5. 

We now proceed to determine the combinations of sfermion masses that would satisfy 
sum rules which are independent of the scalar masses at to, and the D— term contribution. 
We first consider the conventional case. To this end, we note that the right hand side of 
the equations given above that determine the sfermion masses contain a term proportional 
to Yi K, which is proportional to yj(Trym^),and the index i runs over all the sfermions, 
^ = Ql, Ll. ur, (ifi, cr. Wc arc interested in finding those combinations of squared sfermion 
masses where this contribution proportional to K vanishes. Defining this combination as 
Y^- KiUi], where the constants Ki are to be determined, and we have the condition 



J]«,y, = 0. (24) 



We now recall the boundary condition at the GUT scale which reads mf = ml + Yj^gf^D, 
where mo denotes a generic universal soft mass parameter for the sfermion. For the combi- 
nation of squark and slepton masses where no unknown parameters enter, we require that 
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the universal mass term contribution as well as the D term contribution vanish. These 
conditions translate into the equations 

J2Ki^0, (25) 

i 

E '^^^i' = 0- (26) 

i 

We can solve these three constraint equations for the five mass parameters, and in general 
there is a family of solutions. We arbitrarily fix the coefficient of the right handed selectron 
in the sum rules to be —1, and ask for integer values for the parameter Kq for the squark 
doublet. We then arrive at the following two linearly independent sum rules: 

25 

2^1 - ml, - ml = (Ca + 2C2--C^), (27) 
ml + ml - ml -ml = {2C, - ^C,), (28) 
where we have used the notation 

2 ^ t 2 , 2 \ 2 ^ / 2 , 2 \ 

rri/r = -(m,~, + m r ), = -(mx +m,~, )■ 

Q ^2 L '^l' L 2 l ^ 

The essential point is that we have for the conventional embedding Y = X and Y-^ = — Z, 
and have solved the constraint equations (24) and (26). Furthermore, the solutions we 
obtain for flipped embedding is exactly the same. This follows from the fact that the flipped 
embeddding is a special case of the solutions of the constraint equations for the general 
case Y = aX + bZ and Y^ = cX + dZ, where a, b, c, d are arbitrary constants. This is 
not entirely unexpected, since in each of the SU (5) symmetry breaking chains the D-term 
contribution is the same, as each depends only on the gauge contribution. This is one of the 
main results of this paper. It explicitly demonstrates that despite the possibility that the 
model loses its uniqueness in the choice of embedding, the gauge symmetry guarantees that 
certain relations are preserved. 

The two sum rules expressed above, which are valid only when the symmetry breaking 
takes place via SU{5) or flipped SU{5), denote two linearly independent sum rules with 
integer coefficients. The particular choice we have presented is similar to the one in [15]. We 
have chosen to present our sum rules in this manner in order to bring out this feature. The 
reason for our ability to reproduce those sum rules is that the generic non- universality that 
was considered in [15] for the masses of the 10 and 5* has the same structure as is generated 
by D- terms in 5*0(10) unification in the conventional case, where the non- universality is 
now explicitly given by considering the eigenvalues of the corresponding Y-^: 



mlJtc) = rr^Jta) = mUtc) = ml, + gj^D, (29) 

ml, (tc) = ml (tc) = ml, - Sgj.D, (30) 

ml^{tG)=ml,-2glD, (31) 

ml^{tG) = ml, + 2glD, (32) 



at the 5'O(10) unification scale Mq- Here mig and mio are the common soft scalar masses, 
corresponding to the 16 and 10 dimensional representations, respectively of 50(10), at the 
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unification scale. We note here that in the breaking of jS'O(IO) the rank is reduced by one, 
and hence the D-term contribution to the soft masses is expressed by a single parameter D. 
We note from the above that SitQ) = —^gloD. The solution for K is obtained by eliminating 
Ci, C2, C3, m^g, and m^Q from the sfermion mass equations. We get for the conventional 
embedding the result [15]: 

1 10 
X = - (m| - 2mJ, + m\, + m|, -m\ + — sin^ 9wMl cos 2(5) . (33) 



We now consider the eigenvalues corresponding to Y for the flipped embedding. We 
find that the result is identical to that for the conventional embedding, as may be explicitly 
checked. Indeed, the normalization for Y-^ is chosen so as to obtain this result. The conclu- 
sion we draw is that the gauge symmetry of 5'O(10) (i) protects the sum rules irrespective of 
the embedding; (ii) the explicit D— term contributions to the sfermion masses for the two 
embeddings are in fact identical. We note that (ii) implies (i), but not the other way round. 
The stronger result (ii) is of significance to the phenomenology of S0{1Q) models with D— 
term corrections, which is the subject of the next section. 

IV. SPARTICLE SPECTRUM AND CONCLUSIONS 

In order to study the implications of the D- term non-universality for the sparticle 
spectrum, we have carried out numerical integration of the RG equations. We integrate the 
system of coupled renormalization group equations of the gauge couplings, Yukawa couplings, 
the gaugino masses, the scalar mass squared parameters and the soft trilinear couplings 
of the MSSM. We employ well-known results for the unification scale, the value of the 
unified gauge coupling constant and the unified Yukawa coupling for the heaviest generation. 
We present the results of numerical analysis that take into account the effects of the D— 
term non-universality on the boundary conditions for scalar mass squared parameters in the 
conventional and flipped scenarios. For all other parameters we assume universal boundary 
conditions. In other words, starting with values of the common gaugino mass (M1/2), the 
scalar masses at the unification scale (parameterized by mo(= mie = mio) and gf^D), 
trilinear couplings (A), and with the Yukawa couplings having a unified value h at the 
unification scale, and ac, we integrate the set of coupled differential equations down to the 
effective supersymmetry scale of ~ 1 TeV. Of the two minimization conditions for the Higgs 
potential, we recall that one can be written as 

/xf — nl tan^ f3 m% 
tan^ i3 - 1 ~ ~Y 

Proceeding in the well-known fashion [21] of determining tan (3 from the accurately known 
value of the r-lepton mass, and inserting it into the equation above, and using the values of 
the Higgs mass squared parameters determined from the evolution of the parameters yields, 
from the relations /if — "rn^^ + and — + n^, the parameter /i. 

We recall that for the universal case sufficiently large values of the common gaugino mass 
Ml/2 are required to ensure that the gluino is sufficiently heavy, and also fairly large values 
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(34) 



of mo (< M1/2) are required to ensure that the neutrahno is the LSP (in order to prevent the 
hghtest slepton from becoming the LSP). An upper bound on uiq ensues when we require a 
sufficiently large (m^ = /^^ + Keeping these features in mind, we study the effects 
of the D— term on the spectrum. 

In our work we have chosen the sign of D to be positive. This alleviates the problem 
of fine-tuning, inherent in SO (10) unification, by allowing injj^ to evolve to values that are 
negative and larger in magnitude, compared to when the term is absent. (We also note 
here that once solar neutrino parameters are better determined, one would have to include 
the effects of the Yukawa neutrino couphng between Mq and a scale Mr ~ 10^^-^ GeV 
which would also have an effect of alleviating the fine-tuning problem.) We begin with a 
simple illustration for choices of the parameters of the model that would guarantee radiative 
electroweak symmetry breaking, that the lightest neutralino is the LSP and that rriA is 
significantly larger than mz- In Figs.l and 2 we plot the evolution of the five sfermion 
mass squared parameters in a manner which is analogous to [13]. In Fig. 1(a) and 1(b) we 
illustrate the evolution of the mass parameters for the case with D — ior the two lightest 
generations and the heaviest generation, respectively. We have plotted the ratios mj/ml 
vs. t/Log(10), for all the sfermions of interest. We have illustrated this for a typical case of 
Ml/2 = 800, mo = 700, Aq = in units of GeV and the common Yukawa coupling is taken to 
be ht = hb = hr = 2.0. We now add a D— term contribution with gfoD = 300 GeV, which is 
a fairly large contribution. We illustrate this for the lighter generation (Fig. 2(a)) and for the 
heaviest generation (Fig. 2(b)). It is immediately apparent that for the right and left handed 
sleptons of the heaviest generation the renormalization group evolution compensates for the 
initial mass splitting at the GUT scale. As seen clearly in Fig. 2(b), the near degeneracy of 
the two slepton states in the absence of electroweak symmetry breaking makes the mixing 
between them, once SU{2) x C/(l) is broken, significant. Indeed, it is important to observe 
the variation of the mass of the hghest slepton, since it has the tendency to become fighter 
than the hghtest neutralino and to emerge as a candidate for the LSP, and is therefore not 
acceptable. 

We have also plotted in Fig. 3 curves corresponding to important observables. We begin 
by observing the behaviour of /i which is seen to increase significantly as D' = gf^D is 
increased. Note that the parameter // is an important measurable quantity in the chargino 
spectrum. The pseudoscalar Higgs mass m^ which is also plotted in the same figure is 
primarily responsible, apart from tan/?, for determining the properties of the entire Higgs 
spectrum of the model. We have not plotted the value of the mass of the lightest Higgs boson 
since in the regime of uia ^ mz and tan/? ^ 1, its value is a constant of about 120-125 
GeV. It has been pointed out that one can view the rise of uia as a relaxation of fine tuning 
of the parameters in the model [15]. In Fig. 3 the smaller mass eigenvalue resulting from 
the slepton mass squared mixing matrix marked f is plotted. Indeed it rises for a while until 
the mixing effects become pronounced and begins to fall. In particular, the value of D' = 
300 GeV is where the lighter slepton becomes degenerate with the bino whose mass is ~ 350 
GeV when M1/2 — 800 GeV. We have also varied the parameters of the model and find the 
trend to be maintained for reasonable variations. 

To summarize, in this work we have considered what is perhaps the simplest example 
of the D— term non- universality following from 5*0(10) unification. We have obtained sum 
rules for sfermion masses for the conventional and fiipped scenarios in a very general manner 
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and the method may be easily adopted for scenarios based on other groups, e.g., E^, etc. It is 
interesting to note that the same sum rules emerge in conventional and flipped embeddings. 
It turns out that the results are similar to that of non- universality of the type discussed in 
refs. [13,14]. The presence of D— term non- universality assists in rendering the model to be 
technically natural. 
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Figure Captions 



Fig. 1 (a) Evolution of the ratio of tlie sfermion mass squared parameters to the universal 
mass squared parameter, as a function of the momentum scale, for the lightest generations. 
Here we have taken D — 0. The values of the other parameters are given in the text. 

Fig. 1 (b) Evolution of the ratio of the sfermion mass squared parameters to the universal 

mass squared parameter, as a function of the momentum scale, for the heaviest generation. 
Here we have taken D — 0. The values of the other parameters are given in the text. 

Fig. 2 (a) Evolution of the ratio of the sfermion mass squared parameters to the universal 
mass squared parameter, as a function of the momentum scale, for the lightest generations 
with qIqD = 300 GeV. The values of the other parameters are given in the text. 

Fig. 2 (b) Evohition of the ratio of the sfermion mass squared parameters to the universal 
mass squared parameter, as a function of the momentum scale, for the heaviest generation 
with qIqD — 300 GeV. The values of the other parameters are given in the text. 

Fig. 3 Values of mass parameters at the supersymmetry breaking scale, for the choice 
of parameters at the unification scale of M1/2 = 800, mo = 700, Aq — 0, all in GeV, and 
2.0, when D' is varied from to 400 GeV. 
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